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Abstract 

We provide an exhaustive spectral analysis of the two-dimensional peri- 
odic square graph lattice with a magnetic field. We show that the spectrum 
consists of the Dirichlet eigenvalues of the edges and of the preimage of 
the spectrum of a certain discrete operator under the discriminant (Lyapunov 
function) of a suitable Kronig-Penney Hamiltonian. In particular, between 
any two Dirichlet eigenvalues the spectrum is a Cantor set for an irrational 
flux, and is absolutely continuous and has a band structure for a rational flux. 
The Dirichlet eigenvalues can be isolated or embedded, subject to the choice 
of parameters. Conditions for both possibilities are given. We show that 
generically there are infinitely many gaps in the spectrum, and the Bethe- 
Sommerfeld conjecture fails in this case. 



Introduction 

The Hamiltonian H of a two-dimensional Bloch electron in a uniform magnetic 
field B has a highly non-trivial spectral and topological structure depending on the 
ratio of the area a of the elementary cell of the lattice in question and the squared 
magnetic length l 2 M = Hc/\eB\ (here e is the electron charge, c is the light velocity, 
h is the Planck constant). More precisely, denote by 6 the number of the magnetic 
flux quanta through the elementary cell: 6 = c /Ini^. If 8 is a rational number, 
then the spectrum of H has a band structure and for 6 / each vector bundle of 
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the magnetic Bloch functions corresponding to a completely filled Landau levels 
is non-trivial [34] (for a non-zero integer 6 the Chern number of this bundle is 
exactly the values of the quantized Hall conductance in units e 2 /2nh [7,36,40]). 
The most likely conjecture is that the at irrational values of d, the spectrum of 
H has a Cantor structure. This conjecture was recently proved in the case of the 
tight-binding model for the magnetic Bloch electron [5, 6, 37]; in this case H is re- 
duced to the Harper operator in the discrete Hilbert space l 2 (Z). Due to the Lang- 
bein duality [29], the same is true for the Hamiltonian of the nearly-free-electron 
approximation. As a result of this Cantor properties, the diagram describing the 
dependence of the spectrum on the flux 6, so-called Hofstadter butterfly, has a 
remarkable fractal structure, see e.g. [20,22,26]. 

A few is known for the magnetic Schrodinger operator H in the space L 2 (R 2 ), 



where A is the vector potential of B and V is a potential which is periodic with 
respect to the considered lattice. It is proven in this case, that H has a piece of the 
Cantor spectrum near the bottom of the spectrum for a restricted class of potentials 



In this connection, the quantum network models (also called the quantum graph 
models) have attracted considerable interest last time. These models combine 
some essential features of both discrete and continuous models mentioned above. 
On the one hand, the Hamiltonian of a magnetic network model has infinitely 
many magnetic bands of different shape. On the other, the time-independent 
Schrodinger equation for this Hamiltonian can be reduced to a discrete equation. 
S. Alexander was the first who performed this reduction [4] in the framework of 
the percolation approach to the effect of disorder on superconductivity proposed 
by P. G. de Gennes [12]. A very short and elegant derivation of the Schrodinger 
equation for a periodic quantum graph with a uniform magnetic field and a constant 
potential on each edge of the graph is given in [24]. On the mathematical level of 
rigor the relation between solutions of the Schrodinger equation for H on quantum 
graphs r and those for a Jacoby matrix J(H) on the corresponding combinatorial 
graphs was established by P. Exner [15]. Nevertheless, the main theorem from [15] 
is applicable only to the finding of eigenvalues of H distinct from the Dirichlet 
eigenvalues on the edges of T. As to the points of the continuous spectrum, the 
main result of [15] allows an exhaustive analysis only in the case when the direct 
and inverse Schnol-type theorems are known for both H and J(H). 

It is worthy to note that quantum networks are not only a mathematical tool to 
get simplified models of various quantum systems, but in many cases experimental 
devices really have a shape of planar graphs such that the width of the sides is 
much smaller than the parameters of the dimension of length which characterizes 
the quantity in question, e.g., much smaller that the magnetic length, the Fermi 
wave length, the scattering length etc [1,32,33]. In these cases the quantum graph 
models are the most adequate ones for simulating spectral, scattering, and transport 
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V [21]. 
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properties of these devices. 

Here we propose an alternative approach to the spectral analysis of quantum 
graph Hamiltonians based on boundary triples, Dirichlet-to-Neumann maps, and 
the Krein technique of self-adjoint extensions. Such a machinery works effectively 
in many other problems connected with explicitly solvable models [3, 35]. In the 
case of square network lattices with a periodic magnetic field (including a uni- 
form one), an arbitrary L 2 -potential on edges and 5-like boundary conditions at 
the vertices (including the Kirchhoff boundary conditions), we perform an exhaus- 
tive spectral analysis of the network Hamiltonian H. It is proved that the spectrum 
always contains Dirichlet eigenvalues of the edges as infinitely degenerate eigen- 
values of H. The rest part of the spectrum is absolutely continuous and has a band 
structure, if d is a rational number, and is the union of countably many Cantor sets 
placed between Dirichlet eigenvalues, otherwise. Moreover, this part is the preim- 
age of the spectrum of the corresponding lattice Hamiltonian J(H) with respect to 
a many-sheeted analytic function which is just the discriminant of a generalized 
Kronig-Penney operator (i.e. the Sturm-Liouville operator with a Kronig-Penney 
potential). The eigenvalues are isolated if the magnetic flux is non-integer, while 
for integer magnetic fluxes it depends on the electric potential and on the coupling 
at the nodes. 

1 Magnetic Schrodinger operator on the periodic metric 
graph 

Consider a planar square graph lattice whose nodes are the points K m n := (ml,nl), 
(m,n) G I?, where / > is the length of each edge. Two nodes K m n and K p q are 
connected by an edge iff \m — p\ + \n — q\ = 1. We denote the edge between K m n 
and K m+hn by E m ^ r (right), and between K m , n and K mfl+X by E m ^ u (up). Each 
edge E m n r / U will be considered as segment [0, 1] so that is identified in the both 
cases with K m>n , and I is identified with K m+hn for E m ^ r and K m>n+ i for E m ^ u , 
respectively. The phase space of the lattice is 

(m,n)eZ 

The elements of J$? will be denoted as / = (f m ,n,nfm,n,u)> fm,n,r/u £ ^m,n,r/u- 

On 

each edge consider the same electric potential V 6 L 2 [0,/]. 

We assume that the lattice is subjected to an external periodic magnetic field 
orthogonal to the plane, B(x) = (0,0,b(x)) , b G C(M 2 ), b(xi,x 2 ) = b(x\ + l,X2) = 
b(x\ ,X2 + for all x = (x\ ,x%) G M 2 . Denote 

where F is the fundamental domain of the lattice (i.e. the square spanned by £b,o,r 
and Eofl,u) an d represent B(x) = Bq + B p (x), Bq = (0,0,27r<^). The corresponding 
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magnetic vector potential in the symmetric gauge can be written as 

A(xi,x 2 ,x 3 ) = -B Q x (xi,x 2 ,x 3 ) + (Ai(xi,x 2 ),A 2 (xi,x 2 ),0) 

with A ; (xi,x 2 ) = Aj(x\ + /,x 2 ) = Aj{x\ ,x 2 + /), j = 1,2. The presence of the mag- 
netic field leads to non-trivial magnetic potentials on the edges, which are the pro- 
jections of A(x) on the corresponding directions. The magnetic potentials A mnr / U 
on E„ w/u are: 

A m , n , r {t) = (A((ml,nl,0) + (1,0,0)*) , (1,0,0)) = -%nl£, +a r (t), 
A m , n ,u{t) = (A{(ml,nl,0) + (0, 1,0)?) , (0, 1,0)) = nml$+a u (t), 
a r (t) :=Ai(f,0), 
a u {t) :=A 2 (0,0, f€[0,Z]. 
On each of the edges E m n r / U we consider the operator 

-( d V 

with the domain // 2 [0,/]. The direct sum of these operators over all edges is not 
self-adjoint, and in order to obtain a self-adjoint operator on the whole lattice it 
is necessary to introduce boundary conditions at each node. The most general 
boundary conditions involve a number of parameters and can be found, for exam- 
ple, in [25]. We restrict ourselves by considering the so-called magnetic 5-like 
interaction at K m n , 

fm,n,r(ty — ~^ fm,n,u(fi) — fm—\ji,r\f) — ~n fm,n— l t u(f) — • fm,m 
^~dt iAm—l,n,r)fm—l,n,r(fy P (^ 'A m ) «— 1 ,u) fm,n—\,u{f) — (%fm,m 

m,n€Z, 

where a € R, j3 G R \ {0}. These quantities have the following physical meaning. 
The parameter a ^ is the coupling constant of a 5-like potential at each node. 
Introducing the parameter /3 can be treated as considering a more general form of 
the Hamiltonian H: 

H = ^{pj-'?j){*-7*)+y<*y), 

where mjk is the effective mass tensor and /3 is the corresponding anisotropy coef- 
ficient (the ratio of the eigenvalues of the symmetric matrix {m^)). In particular, 
if j8 = 1, one obtains H in the form Q; if in addition a = 0, we get the magnetic 
Kirchhoff coupling. This class of boundary conditions covers main couplings used 
in the physics literature. The self-adjoint operator obtained in this way we denote 
by L. 



4 



Remark 1. It is worthy to note that L is invariant with respect to the so-called 
magnetic translation group Gm [41]. In our case this group is generated by the 
magnetic shift operators z r and % u , 



trfm,n,r/u\t) — ex P 
tufm,n,r/u\t) = ex P 



7ii0[n + 8 u>r / u - 



nid (m + 8 v iJ- 



fm—l,n,r/u\t)) 
fm,n—\,r/u\f)- 



The properties of this group depend drastically on the arithmetic properties of 
6 [10]. In particular, if 6 is irrational, then Gm has only infinite-dimensional irre- 
ducible representations which are trivial on the center of Gm- Therefore, for any 
irrational 6 each point of specL is infinitely degenerate. 



2 Gauge transformations 

To study the spectral properties of L it is useful to use the gauge transformation 

(fm,n,rifm,n,u) — (U mnr (p mnr ,U mnu (p mnu ) given by 

fm,n,r/u(t) exp J A mnr J u ds^j(p mnr j u (t^ . U mnr j u (p mnr j u (t^ 7 

tj w ( t \ — -innl^t+b r (t) ff \ 
u m,n,rYm,n,r\ l ) — c Ym,n,r\ l )i 

tj w ( t \ — Jnml^t+b^t) (t\ 

b r / u (t) ■■= a r/u (s)ds. 

d 2 

There holds U m \ r / u L m ^ r / u U m nj ./ u = ~-^~2 anc ^ ^ e boundary conditions © 
for <p = (<p m ,„ jr , (p m , n ,u) G t/ _1 (domL) take the form 



9m,n = a<Pm,n, (3a) 



where 



(Pm,n • — (Pm,n,r(y) — ^ •Pm,«,w(0) 

= e -^+K, %l _ lnr (l) = L inme+ ^( Pm ^ hu (l), (3b) 

C« :=< ? ) m,^( )+^<«,»( ) 

Kr/u :=b r/u(l)> m,n e Z, 

and 6 := is the number of flux quanta through the elementary cell. Therefore, 
the operator L = U~ l LU acts on each edge as (p m>nir / u i-> — (p£ n r ^ u + V(p m n r j u on 
functions (p satisfying ©, and its spectrum coincides with the spectrum of L. 
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To simplify subsequent calculations we apply another gauge transformation, 
<Pm,n//H = exp (i(mK r + nK u )^ (j> m n r / u . In this way we arrive et an operator acting 
on each edge as <j) m ^ r / u <-> -ty'm,n > r/u +y< t>m,n,r/u on functions = {§ m , n ,r,<$>m,n,u), 
(t>m,n,r/u S H 2 [0,l], satisfying 

<j)' mn = a<j) m ^ n , m,n, e Z, (4a) 
<l>m,n := (j>m,nA°) = \ = e'^^-wCO = ^e iKm6 ^ 1)M (Z), (4b) 

<n:=<n,r(0)+/3< n , M (0)-^^ (4c) 

This operator, which we denote by A, is unitarily equivalent to the initial magnetic 
Hamiltonian. In what follows we work mostly with this new operator. 

At this point we emphasize two important circumstances. First, for usual mag- 
netic Schrodinger operators in L 2 (M 2 ) the spectral analysis for non-zero but peri- 
odic magnetic vector potentials (i.e. with the zero flux per cell) essentially differs 
from that for the Schrodinger operators without magnetic field. Even the proof of 
the absolute continuity of the spectrum is non-trivial [8,39]. In our case, the oper- 
ator on the graph with a periodic magnetic vector potential appears to be unitarily 
equivalent to the operator without magnetic field. Second, for the usual magnetic 
Schrodinger operators the bottom of the spectrum grows infinitely as the flux be- 
comes infinitely large. In our situation, the spectrum is 1 -periodic with respect to 
the magnetic flux 8, as changing 6 by 6 + 1 in © obviously can compensated by 
a unitary transformation. Such periodicity leads to the so-called Aharonov-Bohm 
oscillations in the corresponding physical quantities. 

Proposition 2. The operator A is semibounded below. 

Proof. Let <p £ domA. Using the integration by parts and changing suitably the 
summation order one obtains 

<0,A0) = £ ((0, w , -^m,n,r + V^m,n,r) + {<l>m,n,u,-<l>m,n,u+ V( l>m,n, U )) 
m,n 

= E(^/(0)^ )r (0)-^, r (/)^^(/) + ^, il (0)^^(0)-^ ll (/)^, )fl)1| (/) 
m,n 

+ J q (\ K,,n,r 1 2 + V\ <j> m ,n,r \ 2 )dx + J^(\ §' m n u | 2 + V | m> „,„ | 2 ) dx) 

= L{ [ {W m .n,r\ 2 + V\K,nA 2 )dx+ f fl^J 2 + V\^ )U \ % )dx 
tn n 

+^(< n/ (0) + j3<„, M (0)-e^^^ 

= L { J Q ( I ^'m,n,r ? + V \ <l>m,n,r | 2 ) dx + jf ( | ^ J 2 + V \ (j> m , n ,u \ 2 )dx + ^ n ^ m ^ 

= l{ f {W n ,n,r?+y\^A 2 )dx+ /' ( + V\ $m,n,« \ 2 )dx+ 0f| (j> m , n | 2 } . 
m,n •'0 •'0 
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Now choose c € (0, 1) and CgM with 

\a\\h(0)\ 2 < / ( c\ti ) 2 + (V + C)\h\ 2 ) dx for all h £ H l [0,1] 
Jo 

(the existence of such constants follows from the Sobolev inequality), then 
|«||<M 2 = Ml4w(0)| 2 < J o {cW m ^ r \ 2 + {V + C)\^ n A 2 )dx, 



and 

(0,(A + C)0)>£ /'{(l-c)|<„ )r | 2 + |^J 2 +(y + C)|^„, M | 2 }jx>0. □ 

3 Boundary triples 

Here we describe briefly the technique of abstract self-adjoint boundary value prob- 
lems with the help of boundary triples. For more detailed discussion we refer 
to [13]. 

Let S be a closed linear operator in a Hilbert space Jif with the domain domS. 
Assume that there exists an auxiliary Hilbert space & and two linear maps F, V : 
domS — ► & such that 

• for any f,g e dom <> there holds (f,Sg) - (Sf,g) = (rf,Pg) - (Pf,r g ), 

• the map (r, H) : dom S —>■ W © & is surjective, 

• the set ker T n ker V is dense in 

The triple (^,r,r') with the above properties is called a boundary triple for S. 

Example 3. Let us describe one important example of boundary triple. Let V G 
L 2 [0,Z] be a real-valued function. In = L 2 [0,l] consider the operator 

d 2 

S = -— + V, domS = H 2 [0,l], (5) 



then one can set 



If an operator S has a boundary triple, then it has self-adjoint restrictions (see 
proof of theorem 3.1.6 in [19]), hence 5* is a symmetric operator. For example, if 
T is a self-adjoint operator in W, then the restriction of S to elements / satisfying 
abstract boundary conditions F'f = TTf is a self-adjoint operator in J^, which we 
denote by Ht- Another example is the operator H corresponding to the boundary 
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conditions Tf = 0. One can relate the resolvents of H and Hj as well as their 
spectral properties by the Krein resolvent formula, which is our most important 
tool in this paper. 

Let z ^ spec//. For g G 5f denote by y(z)g the unique solution to the abstract 
boundary value problem (S — z)f = with Tf = g (the solution exists due to the 
above conditions for T and T'). Clearly, y(z) is a linear map from 'S to M '. Denote 
also by Q(z) the operator on given by Q(z)g = T'j(z)g; this map is called the 
Krein function. The operator- valued functions 7 and Q are analytic outside spec//. 
Moreover, Q(z) is self-adjoint for real z. 

Proposition 4. (A) (Proposition 2 in [13]) Forz $ spec //U spec Hj the operator 
Q(z.) — T acting on <S has a bounded inverse defined everywhere and the 
Krein resolvent formula holds, 

(H-z)- l -(H T -z)- l = Y(z) (g(z) - T) ~ V (z). 

(B) The set specHj \ spec// consists exactly of real numbers z. such that G 
spec(<2(z)-r). 

( C) (Theorem 1 in [17]) Let z, G spec//^ \ spec//, then z is an eigenvalue of Hj 
if and only if is an eigenvalue of Q{z) — T, and in this case y{z) is an 
isomorphism of the corresponding eigensubspaces. 

This statement is especially useful if the spectrum of H is a discrete set and the 
spectrum of Hj is expected to have a positive measure, because one can describe 
the most part of the spectrum of Hj in terms of Q(z) — T. 

Example 5. Consider the example given by (|5j and ©. The corresponding Krein 
function s(x) can be obtained as follows. The restriction D of S given by Tf = is 

Df=-f" + Vf, domD = {/G// 2 [0,Z]:/(0)= /(/) = 0} (7) 
In what follows we denote the eigenvalues of D by jl^, k = 0, 1,2. . . , /Xq < Ml < 

M2 < 

Let two functions U\,U2 G H 2 [0,l] satisfy 

cv ,„ , «i(0;z)=0, ui(0;z) = l, 
wi,w 2 Gker(S-z), (8) 
"2(0;z) = l, m 2 ( ;z) =0 - 

Clearly, for their Wronskian one has w(z) = u' 1 (x;z)u2(x;z) — u\{x;z)u' 2 {x;z) = 1 
Both u\, U2 as well as their derivatives with respect to x are entire functions of z. 
Let z $ spec D, then any function / with —f" + Vf = zf can be written as 

= 77-^ ui(x;z)+f(0)u2(x;z), 

u\{l\z) 
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and the calculation of f'(0) and —/'(/) gives 



-U2(l;z 
w(z 



-U2(l;z) 
l 



-«i(/;z), 



(9) 



It can be directly seen that s(z) is real and self-adjoint for real z. Clearly, the matrix 
s has simple poles at jU^, which are at the same time simple zeros of u\ (Z;z). More 
precisely, by the well-known arguments, see e.g. Eq. (1.4.13) in [30], there holds 



du\(l;z) 



dz 



Ui 



(l;Hk) [ u\{s, 
Jo 



and U2(l',Hk) 7^ due to i^i(l;Hk)u2(l',Hk) = w(fik 



4 Reduction to a discrete problem on the lattice 

To describe the spectrum of A we use the Krein resolvent formula. Denote by IT the 
operator acting on each edge as tp mjur / u i— > —<j>^ n r t u + Vty nhn r / U on functions satis- 
fying only the condition d4bb . Clearly, for such functions the expression ty' m n given 
by d4cl makes sense. This operator is not symmetric, as it is a proper extension of 
the self-adjoint operator A. 

Proposition 6. The operator IT is closed and the triple (/ 2 (Z 2 ),r,r'), 

T : domll 9 = {§ m , nj r,tym,n,u) ^ (0m,n) S f(Z 2 ), 

T : domn 9 = (<t> mjl , r ,<t> mjl ,u) » ««) G l 2 (Z 2 ), 
is a boundary triple for IT. 

d 2 

Proof. Denote by S the direct sum of operators — + V with the domain H [0, /] 
over all edges E mn r / U . Clearly, S is a closed operator, and the functionals 

gm,n,l(4>) = 0m,n,r(O) 
gm,n,2(<l>) = 0m,n,r(O) 
gm,n?{§) = 0m,n,r(O) 

are continuous with respect to the graph norm of S. Therefore, the restriction of 
E to functions on which all these functionals vanish is a closed operator, and this 
is exactly IT. For any G domll the inclusions r0,P0 G l 2 (I?) follow from the 
Sobolev inequality, and both T, F' are continuous with respect to the graph norm 
of n. 



(0), 

-e- ime <j) m -u,Al), 
-\e inmQ $ m , n -,, u (l) 
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Let 0, y G domlT, then the integration by parts gives 

(4>,nv)-(rty,v) 

= £ ((<!>m,n,i,-Vm,nj) ~ {~K,n,v ¥m,n,i)) 
m,n^Z,i=r,u 

= L {i,»,r(0)i/ ra , v (0) + i^(0)/3<„, M (0) 

-^m-l,n,r(/)v4-l,n,r(0-^m,n-l,«(/)/3v4,n-l,«W 
" C^(0)Wn,n,r(0) " ^ (0) 

+ C-l,»,r( Z )W»-l.n,»-(0 + £ (O0Vm,«-l,«(O} 
= X {^<.«,r(0)+/3^<,„ jM (0) 

-C,n,r(0) W" - j8 (0) Vm,n 

+ Vk» + PK^)e- inme ¥m,n } 

= £ {^«„,r(0)+j3<„, M (0) 

-e- im W m -^l)-^ m W m ,n-i,u{l)) ~ 

-e- ,m9 €-UnA l )-P ei * md €,n-a i )) Ym,n} 

Now we verify the surjectivity condition. Choose functions fo,fi G H 2 [0,l], with 
/o(0) = /{(0) = 1, f (0) = /i(0) = 0, /^(Z) = 0, j,k = 0, 1. Let g,g> e Z 2 (Z 2 ). 
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For any p,q G Z denote 



h p , q ,i (x) = g p , q fo(x) + -^/i (x), 



h p , q , 3 (x) =g P , q e inqB f (l-x)+e- i ^ e6 -^f l {l 
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h p , qA (x) = (3g p , q e-^ d Ml - x) + ^-fi(l-x), 
then /i p , 9 ,y G // 2 [0, /] , j = 1 , . . . , 4, and these functions satisfy 

W(O) = 4 W(0) =e- iKqe h p ^{l) = ^ np6 h MA {l) = g M , 

ti p , q ,i(0)=Ph' p>q!2 (0) = -e^h'^l) = -^e-^h' MA (l) = 8 -f, 

h p , q ,i(l) = h p ^ 2 {l) = h p ^{0) =h PtqA (0) 
= h' p , q ,i (I) = h' m2 {l) = h' pjq>3 (0) = ti MA (0) = 0. 

Define = G with 

=h p , q ,i, = h p ^ 2 , = h p>qj3 , $ { p P q q \ ;U = h MA , 

§mn] = f° r a ^ otrier "!,ngZ and i = r,u. 

Clearly, by construction 0^ G doml! and there holds {T^ p ' q) ) mn = ty'tfjfi = 

g P , q 8mp8„ q and {^'^> {p ' q) ) mn = {<t> ip ' q) )' m ^ = g'p, q SmpS nq , m,neZ. It is easy to see 
that the series = £ mn ^ m,n ' converges in the graph norm of IT, hence G domll. 
Since // 2 [0,Z] is continuously imbedded in C 1 [0,Z], we have F(j) = Y,mn^^ m = g 
and T'(j) = Em«r / ^ = <?'■ The surjectivity condition is proved. 

It remains to note that the set ker T n ker V contains the direct sum of Cq(0, I) 
over all edges and is obviously dense in J$? . □ 

The operator A is the restriction of IT to the set of function <p satisfying Vtp = 
aTip. Consider another self-adjoint extension Ilo given by Tip = 0. Clearly, Ilo is 
exactly the direct sum of the operators D from (Q over all the edges E m n r / U . In 
particular, spec Ilo = spec D. 

Let z ^ specD, g G l 2 (I?) and \j/ g be solution of (II — z)y g = satisfying the 
boundary condition r\j/ g = g. Consider the corresponding Krein function Q(z) : 
I 2 (Z 2 ) 3 g i— > r'y/'g G Z 2 (Z 2 ). Application of proposition 0] gives the following 
implicit description of the spectrum of A. 

Proposition 7. A number z G R \ spec Ilo = R \ specD lies in spec A iff G 
spec [Q(z)-a]. Such Z is an eigenvalue of A iffO is an eigenvalue of Q{z) — a. 
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Therefore, outside the discrete set spec ITo = spec D we can reduce the spectral 
problem for A to the spectral problem for Q(z) — a. Let us calculate Q(z) more 
explicitly; actually this is our key construction. 



(10) 



Proposition 8. For z ^ spec D there holds 

Q(z) = (1 +/3 2 ) {su(z) + s 22 (z)) + Sn (z)M(9,P). 
where M (6, (5) is the discrete magnetic Laplacian in Z 2 (Z 2 ), 
(M(9,P)g) mn = e im9 gm+l , n + e- ime g m - lM +l5 2 (e^ 9 g m , n+1 + e i70nd g m ^) , 

Proof. Note that for = {§ m ^ n ,r, <l>m,n,u) i n tne notation of proposition [6]there holds 



/ f<l>m,n,M\ \ 
0m,«,r(O 
0m,«,«(O) 

\ \<t>m,n,u( l ) / / ( m „„) G Z 2 



with operators B : I (Z ) 



/ / ^,„,r(0) \ \ 

K,n,u(°) 

(m.n) 



B : 



( \ 

h (2) 

h {3) 
,l m,n 



(h%n + e- i ™%^ n + f5h% n + l3e i ™%Z- 1 ), 



Z 2 (Z 2 ) andC: / 



(3) 



Z 2 (Z 2 )®C 4 given by 



Jnm6 h (4) 



and 



C ■ (gm.n) 



Jnn9 _ 
e gm+l,n 



figm,n 



—iTlmQ , 



i,n+l/ J 



Cf. ©. 

Let g 6 Z 2 (Z 2 ). For z ^ spec S, finding the solution <p with (H — z)<p =0 and 
r0 = g reduces to a series of boundary values problems for components of 0, 



f d £ \ 
[--^2+V(t)-z)^ m , n , r/u (t)=0, 

Cg, 



( /4w(o)\ \ 

0m,n,«(O) 
V \ <Pm,nA l )J J 
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and 







fsu(z) 







\ 


1 (j>m,n,r(°)\ 


-0m,n,r(O 




521 (z) 


522 (z) 










C,n,„(0) 










s n(z) 


*12(z) 


0m,n,«(O) 


\-<S>'m,n,u(l)) 









521 (z) 


s 2l{z)) 


\§m,n,u{l)) 



For <2(z)g = r'0 one has 



r> = 5 



/ / C,„,,(o) \ \ 

-0m,n.r(O 

vv-<„,„(o/y 



Therefore, <2(z) = BK(z)C, where A'(z) is a linear operator on / 2 (Z 2 
matrix 

/All(z) (z) 

521 (z) 5 22 (z) 



VV o o 

In other words, for any g € Z 2 (Z 2 ) one has 



K(z) = diag 



\\ 


511 (z) 512 (z) 
521 (z) S 22 (z)J J 



I ( 



Cg 



Jnn8 „ 

« gm+l,n 



K(z)Cg 



Pgm.n 

V \Pe-^ me g m , n+l J ) 

(( Su(z)gm,n + e inne s n (z)gm+l,n \\ 
S2l{z)gm,n+e mn6 S22{z)g m +l,n 

P{sn(z)g m ,„ + e- i7lm0 s l2 (z)g m ,n + i) 

v Vi3(52i(z)^m,„+e-'' 7rme s22(z)^ ! ,, 1+ i)y y 



4 with the 



and, finally, 

= (1 +/3 2 ) (511(Z) +5 2 2(z))g m ,„ + ^' ie 5l2(z)gm + M+e- , ' Jrne 52l(z)^-l > » 

+ j 8V^ e , 12 (^ m ,„ +1 + j3V D "%(z)fc B , ll _i. (11) 
As can be seen from ©, there holds si2(z) = 521 (z) an d becomes exactly 

<EoJ. □ 

Corollary 9. A number z € M \ specD Ziej in f/ze spectrum ofL iff 

e spec [(1 +/3 2 ) (s n (z) +522 (z)) - a + 5 12 (z)M(0,/3)] . 
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5 Spectral analysis 

To describe the spectrum of A we need some additional information on the Krein 
matrix s(z) from ©. 

Proposition 10. The matrix s(z) has the following properties: 

(A) s l2 (z) ^0 for all spec D. 

(B) For any a E R the function 

sn(z) s n (z) 
can be extended to an entire function. 

(C) The function y-— ^ T](z) is the discriminant of the (generalized) Sturm- 



Liouville operator 



d 2 

-—+W(t)+W KP (t), (13) 



where W is the periodic extension ofV, W(t + nl) = V(t), t E [0,1), n E Z, 

(X 

and Wkp is the Kronig-Penney potential Wrp(0 = ^ 8(t — kl); such 

operators P also are called Kronig-Penney Hamiltonians. 

(D) There holds T](}X k ) < -2(1 + j8 2 ) for even k and T}(}l k ) > 2(1 + j8 2 ) for odd 
k. (Recall that \x k are the eigenvalues ofD.) 

(E) For all real z with |tj(z)| < 2(1 + j8 2 ) there holds T]'(z) ^ 0. The function 
f] has no local minima with f] = 2(1 + /3 2 ) and no local maxima with f] = 

-2(1 +n 

Proof. Recall that s(z) is given by © with u\, ui from (J8J). There holds s\i(z) = 

— - — r- and sn(z) ^ for all z ^ specD since z *-* «i(/;z) is an entire function. 
u\(l;z) 

This proves (A). 

Substituting Q for z ^ specD in (fT2b one arrives at 

tj ( Z ) = (1 + j8 2 ) («; (Z;z) + «a(Z;z)) + a«i (Z;z) 

and T] obviously has analytic extension to all points of specD. This proves (B). 

To understand the meaning of v\ look at the operator (fT3l . This operator acts 
as / i-» -/" + Wf on functions f£H 2 (R\ IZ) satisfying 
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Let y\, y>2 be two solutions of (P — z,)y = with y\ (0+;z) = y' 2 (0+;z) = and 
y' l (0+;z) =y2(0+;z) = 1. Consider the matrix 



M(z) 



y[(l+;z) y' z (l+;z) 
yi{l+;z) y 2 (i+;z) 



It is well-known that the spectrum of P consists exactly of real z satisfying rrM(z) = 
y[{l+;z) +y2(l+',z) G [-2,2], see e.g. [16,23]. The function trM(z) is called the 
discriminant or the Lyapunov function of P and plays an important role in the 
study of second order differential operators; if a = 0, the study of this function is 
a classical topic of the theory of ordinary differential equations, see e.g. [11,30]. 

On the other hand, note that on the interval (0,/) the solutions y\ and y 2 coin- 
cide with u\ and u 2 from ©, respectively. In particular, yij,(l—',z) = "1,2(^2) and 
y\ 2 {l—',z) = u'j 2(^2) Therefore, taking into account the boundary conditions (fPfl) 
we can write M(z) in the form 

/ OC OC \ 

u\{l;z) + u\(l;z) - — ■^ 2 u 2 {l;z) + u' 2 (l;z)\ 



M[z)= 1+/3 2 1V ' J 1V ' ; l+p 2 , . 

\ u\{l;z) u 2 (l\z) J 

and trM(z) = - ^_^ 2 ' { ]{z), which proves (C). The items (D) and (E) describe typi- 
cal properties of the discriminants of one-dimensional periodic operators. 

To prove (D) note that for any k one has u\{l,pLu) = 0, and T7 (jU^) = (1 + 

fi 2 )((u\{l;Hk) + U2(l\jJLkj), i-e. - ^^ ^iHk) coincides with the value of the dis- 
criminant of the classical periodic Sturm-Liouville problem (a = 0), for which the 
requested inequalities are well known, see e.g. Lemma VIII.3.1 in [11]. 

The first part of (E) is known for much more general potentials, see e.g. 
Lemma 5.2 in [23]. As for the second part, local maxima with v\ = —2(1 + j8 ) 
and local minima with r\ = 2(1 + j8 2 ) would be isolated eigenvalues of P, which is 
impossible, because the spectrum of P is absolutely continuous [31]. □ 

Therefore, up to the discrete set specD the spectrum of A is the preimage of 
specM(0,/3) under the entire function 17. The operator M(0,j3) is very sensi- 
tive to the arithmetic properties of 6 and is closely related to the Harper operator, 
cf. [38]. The nature of the spectrum of M(8, j8) in its dependence on 6 is described 
in the following proposition, which summarizes theorem 2.7 in [38] (item A), the- 
orem 4.2 in [38], theorem 1.6 in [6] and main theorem in [5] (item B), and theorem 
2.1 in [9] (item C). 

Proposition 11. (A) The operator M (6 , j8) has no eigenvalues for all 6 and /3. 

(B) If 6 is irrational, the spectrum ofM(6,[5) is a Cantor set. If in addition, 
j8 = 1, the spectrum has zero Lebesgue measure. 

(C) For non-integer 6 there holds ||M(0,j8)|| < 2(1 + j8 2 ). 
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The previous discussion gives a description of the spectrum of A in E \ spec D. 
Let us include spec D into consideration. 

Proposition 12. There holds specD C spec A. Moreover, each fi^ G specD ij a« 
infinitely degenerate eigenvalue of A. 

Proof. Consider an eigenvalue of D and the corresponding eigenfunction / with 
/(0) = land let a :=/(/). 

Let be rational. Take MeZ such that BM G 2Z. Let p,q £Z. Denote by 
the function from whose only non-zero components are 

<PpM+j,qM,r = fiO J f, <PpM,qM+j,u = —G J f, 

7 = 0,...,M-1. 

Clearly, G domA and — ^ n , + (V — ^k)^m,n,r/u = for all m,n £ Z. There- 
fore, is an eigenfunction of A with the eigenvalue /i^. As p and q are arbitrary, 
one can construct infinitely many eigenfunctions with non-intersecting supports. 
Therefore, each is infinitely degenerate in spec A. 

Now let 6 be irrational. We use arguments similar to the Schnol-type the- 
orems [28]. For each «6Z put 0_„ jn> , = pe i7ln9 f and 0_„ AM = -e i7ln9 f. The 
"chain" constructed from these components does not belong to J^, but satisfies 
the boundary conditions ©. Moreover, — r , u + (V — pLk)ty n -n,r/u = f° r au 
«, i.e. this chain is a "generalized eigenfunction" of A. 

Take (p G C°°[0,l] with <p(0) = <p'(0) = and <p(Z) = p'(Z) = 1. For any N G N 
construct i^ (A,) G Jt? such that 

V^w/« = 4>-n,n,r/u if H < tf, 

(/V) . (iV) . 

¥-N,N,r — 9<P-N^,r, ¥n,-N,u ~ 

vff2 j = f° r ai l other m,n£Z and i = r,u. 

Clearly, y/"^ G domA for any AT and > v 7 ^!!/!!- Moreover, the only two 

non-zero components of g^ N > = (A— are 



-N,N,r 



= - j 8^(«p" / + 2<p7 / ) 1 

and 

= ( " <P"f ~ 2 <P'f ~ <Pf + ( V - «*)*>/) = Wf + 2 <P'f) ■ 



Therefore, \\gW\\ = \\(A- ^ k )\j/^\\ = y/l + p 2 \\(p"f + 2(p'f\\ = C and 

r ||(A-^)y W H ^ r C 
hm - — 7— 11 < hm 



at-h- \\V {N) \\ V2N\\f\\ 
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which means that G spec A. Let us show that is an eigenvalue of A. By 
proposition [TTlC) one has ||M(0,/$)|| < 2(1 + /3 2 ). Recall that the spectrum of A 
outside specD is the preimage of specM(0,/3) under the function v\ and, due to 
proposition HU D), does not contain /i^. As is an isolated point of the spectrum, 
it is an eigenvalue of A, which is infinitely degenerate according to the arguments 
given in remark ^ D 

Now we state the main result of the paper. 

Theorem 13. The spectrum of A is the union of two sets, 

specA = £ U£, £ = specD, 1 = T7 _1 (specM(0,j3)), 

and has the following properties: 

(A) The discrete spectrum is empty and the point spectrum coincides with £o 

(B ) The set £ is non-empty, moreover, the intersection [jU* , jJ-k+ 1 ] H£ is non-empty 
for any lc 

( C) For rational the singular continuous spectrum of A is empty and the abso- 
lutely continuous spectrum coincides with £ and has a band structure. 

(D) For irrational 0, the spectrum of A is infinitely degenerate. The part £ is 
a closed nowhere dense set without isolated points, and £n (jUfc, Jlfc+i) is a 

Cantor set for any k = 0,\,2, If additionally j8 = 1, then the spectrum of 

A has no absolutely continuous part and the singular continuous spectrum 
coincides with £ 

Proof. Proposition 1721 shows that £o C spec A. The spectrum of A outside specD 
is described by corollary [9] and, in virtue of proposition [ToF B ) . coincides with £. 

(A) Propositions [8j [TOl B). and ITiT A) show the absence of eigenvalues of 
Q(z) — CC. In virtue of proposition the operator A has no point spectrum in 
M\specD for any 6. Therefore, due to proposition [T2l the point spectrum coin- 
cides with spec D, and all eigenvalues have infinite multiplicity. 

(B) The trivial estimate ||M(0,/$)|| < 2(1 + j8 2 ) implies the inclusion 
specM(0,j8) C [-2(l+j3 2 ),2(l+j3 2 )]. The assertion follows now from proposi- 
tionUOtD). 

(C) Let be rational. Take N G Z with A^0 G 2Z. The operator A appears to 
be invariant under the shifts (tj) m ^ r , $ m ,n,u) ^ ($m+kN,n+lN,r,$m+W,n+lN,u), k,leZ, 
i.e. is Z 2 -periodic and, therefore, the absence of singular spectrum for A follows 
from the standard arguments of the Bloch theory, see e.g. theorem 11 in [28]. 
Therefore by (A) £ coincides with the absolutely continuous spectrum. The spec- 
trum of M(0,j8) consists of finitely many bands, so is Tj _1 (specM(0,j3)) between 
any two Dirichlet eigenvalues. 

(D) Now let be irrational. The infinite degeneracy of spec A follows from 
the arguments of remarkfj The set of z G M for which |t](z)| < 2(1 + j8 2 ) consists 
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of non-intersecting finite intervals /„. Moreover, due to proposition [TOlD ) there 
is exactly one such interval between any two subsequent eigenvalues of D. Put 
J n = I n . Note that U7„ contains all points z with <2(l + /3 2 ). Due to propo- 

sition QUE), the restriction of v\ to J n is a homeomorphism of J n on the segment 
[-2(1 + j3 2 ),2(l + j8 2 )]. Therefore, the preimage 

Kn := (T]| 7 „)" 1 (specM(0, j 8)) cJ n 

is a Cantor set as is true of specM(0,/3). Moreover, the intersection of any two 
of sets K m is empty, which follows from proposition [TlT C) and proposition [TOlD ). 
Therefore, the set UK n , which coincides with £, is also closed, nowhere dense, and 
without isolated points. 

If p = 1, then the sets K n are of zero Lebesgue measure due to proposi- 
tion ^nB). Since (t]|/„) 1 are real analytic, the set £ = UK n is also of zero 
Lebesgue measure. Such a set cannot support absolutely continuous spectrum and 
does not intersect the point spectrum due to (A), therefore, £ is the singular con- 
tinuous spectrum. □ 

We formulate several corollaries in order to answer the following natural ques- 
tions: 

• Are the eigenvalues of A isolated or embedded in the continuous spectrum? 

• Is the number of gaps in the spectrum finite or infinite? Note that the rank 
of the lattice defining the magnetic translation group is equal to 2, therefore, 
one can expect the validity of the Bethe-Sommerfeld conjecture for 6 = 0. 

We emphasize that these questions are rather non-trivial even for lattices with- 
out any potentials; for example, rectangular lattices with 5-boundary conditions at 
the nodes can have very different properties depending on the coupling constants 
and the ratio between the edge lengths [14]. We will see that the introduction of 
scalar potentials on edges provides a mechanism of gap creation similar to the so- 
called decoration [2]. 

Let us consider first the case without magnetic field in greater details. 

Corollary 14. Let 6 be integer. 

(A) The part £ of the spectrum of A. coincides with the spectrum of the Kronig- 
Penney Hamiltonian P from dl3t . In particular, if there are inifinitely many 
gaps in the spectrum ofP, then A has the same property. 

(B) IfV is a convex smooth function whose derivative does not vanish, then all 
gaps are open for any a, the spectrum ofP and hence also of A has infinitely 
many gaps, and all }Xk are isolated in spec A. 
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(C) Let the gap of P near fa be closed for a = 0. Then fa is an embedded 
eigenvalue of A. for all a. In particular, fa lies on a band edge for a 7^ 0. 
For V = and a =0 all gaps are closed and all fa are embedded into the 
continuous spectrum. 

Proof. (A) In this case one has specM(0,jS) = [-2(1 +/3 2 ),2(1 +/3 2 )] and the 
set £= 77 _1 2(1 +j6 2 ),2(l +/3 2 )]^ coincides with the spectrum of Pby propo- 
sitionEIC). 

(B) Denote v(z) = u\(l;z) + U2{l',z), where u\ and U2 are the special solutions 
from example |5] Clearly, v is the discriminant of the periodic Sturm-Liouville 

d 2 

operator Q := T + W with W from Jl3l . If a = 0, then P = £2 and n(z) = 

(1 + j8 2 )v(z) for all z. Let V be smooth convex with V' / and a = 0, then it 
is proved in [18] (see lemma 3 and theorem 2 therein) that all gaps of P are open 
and that fa do not belong to spec/ 3 = E, which means that |v (/!&)! > 2. One has 
77 (fa) = (1 + fi 2 )v(fa) + au\(l\fa) = (1 + p 2 )v(fa), and the gap remains open 
for all a ^0 as \r)(fa)\ > 2(1 + j3 2 ). 

(C) The case with V = and a = is obvious. If the gap near fa is closed, 
then v(jU)t) = ±2, T7(^) = ±2(1 + j8 2 ), and ju* G I. Moreover, v'(jUjt) = 0. As 
d z u\(l;fa) 7^ (see example^, one has T\'(fa) / for a / 0, which means that 
T7=f2(1+/3 2 ) changes the sign at fa. This means that there is a gap near fa. □ 

The case of a non-trivial magnetic field can be treated by a simple norm esti- 
mate. 

Corollary 15. If 6 is non-integer, then the spectrum of A has infinitely many gaps, 
and all fa lie inside the gap. 

Proof. In this case the set £ does not contain fa due to propositions [TOl D) 
andEllD). □ 

6 Concluding remarks 

One can easily modify the above proof to many other boundary conditions at the 
nodes K„ un , for example, for the so-called ^'-boundary conditions, 



(Pj = a fm,m m,n£Z, 

A suitable boundary triple would be (/ 2 (Z 2 ),ro,rQ) with To = —I* and r' = T 
with r, r' from proposition |6] Instead of the operator D from Q one should deal 
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with the Neumann realization, Nf = -f" + Vf, domN = {/ e H 2 [0,l],f'(0) = 
f'{l) = 0}, and the corresponding constructions can be repeated almost literally. 

The approach can be extended to the analysis of more general periodic mag- 
netic systems, for example, for more complicated combinatorial structures, for 
nodes and edges with geometric defects or measure potentials, or with the spin- 
orbital coupling taken into account. Such systems will be considered in details in 
forthcoming publications. 
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